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Abstract: The Fermi surface calculated within the rotating antiferromagentism theory 
^ 2 undergoes a topological change when doping changes from p-type to n-type, in quali- 

^ s tative agreement with experimental data for n-type cuprate Nd2_a;Cea;Cu04 and p-type 

^ 4 La2-a;Sra;Cu04. Also, the reconstruction of the Fermi surface observed experimentally 

^ 5 close to optimal doing in p-type cuprates, and slightly higher than optimal doping in the 

Q ? overdoped regime for this n-type high-T^ cuprate is well accounted for in this theory, 

, '^, 7 and is a consequence of quantum criticality caused by the disappearance of rotating 

g antiferromagnetism. The present results are in qualitative agreement with the recently 

9 observed quantum oscillations in some high-T^ cuprates regarding the change in the size 

of the Fermi surface as doping evolves and the location of its reconstruction. This paper 
presents new results about the application of the rotating antiferromagnetism theory to the 
psj^? study of electronic structure for n-type materials. 
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>i4 breaking; Fermi surface reconstruction 



15 1. Introduction 



16 The topology and doping dependence of the Fermi surface (FS) of high-temperature superconductors 

17 (HTSC) are currently highly debated. One of the reasons for this debate is due to some contradicting 

18 data collected from angle-resolved-phtoemission spectroscopy (ARPES) experiments, which do not 

19 seem to see any FS reconstraction due to the occurrence of a symmetry breaking, and data collected 

20 from magnetoresistance measurements, characterized by quantum oscillation like Shubnikov-de Haas 

21 (SdH) Oscillations, which reveal that the FS undergoes a topology change due to some sort of symmetry 
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22 breaking. Since no long range order has been observed so far in underdoped HTSCs we proposed earlier 

23 that the FS reconstruction is caused by the hidden rotating antiferromagnetic order. 

24 In the present work we support this proposal by new results for n-type doping case and argue in favor 

25 of the FS of HTSCs undergoing topology reconstruction at specific doping levels in the framework of 

26 rotating antiferromagnetism theory (RAFT) [1]. We compare the evolution of the FS with doping in 

27 p-type and n-type HTSCs obtained in this theory, and discuss it in connection mainly with available 

28 experimental data for n-type material Nd2_2:Cea;Cu04 and p-type one La2_aSraCu04. It is found that 

29 the change in the topology of the FS as one goes from the p-type cuprate to n-type material is well 

30 accounted for in RAFT. In the low-doping limit (underdoped regime) RAFT yields a small almost 

31 square FS centered around (vr, 0) points for n-type Nd2-a:Ce2;Cu04 in qualitative agreement with SdH 

32 oscillations results regarding the existence of a FS in the form of small pockets [2]. A careful look 

33 at the data of Armitage et al. [3] in Fig. 3 of their work reveals a trend qualitatively consistent with 

34 our findings for n-type material Nd2_a;Cea:Cu04 regarding the evolution of spectral weight away from 

35 (tt, 0) to form a larger FS as doping increases. A FS in the form of stretched elliptic pockets nearby the 

36 (vr/2, 7r/2) points is however likely for p-type La2_a:Sra,Cu04. Indeed, Fig. 5 of the ARPES work by 

37 Yoshida et al. [4] shows nicely the evolution of the FS with doping from what we interpret as stretched 

38 small pockets in the underdoped regime to large countours in the overdoped regime. RAFT reproduces 

39 qualitatively well the FS evolution with doping for this p-type material. Note that p-type cuprates were 

40 examined using RAFT in Ref. [5]. However this is the first work based on RAFT, which deals with 

41 the electronic structure in an n-type cuprate. In RAFT, for both p-type and n-type materials, the critical 

42 value of doping where FS reconstruction occurs is given by the value where rotating antiferromagnetism 

43 vanishes. In p-type materials, this value coincides practically with optimal doping, but in n-type case, it 

44 occurs in the overdoped regime beyond optimal doping for superconductivity (SC). 

45 This paper is organized as follows. In Section 2, RAFT is extensively reviewed. In Section 3 the 

46 rotating antiferromagnetic and superconducting parameters are calculated as a function of doping and 

47 temperature. In Section 4, the doping dependence of the electronic stucture is calculated and compared 

48 to experimental data. Energy spectra versus wavevector are calculated for several doping levels, and the 

49 FS is calculated in this section using the occupation probability for doping levels in n-type and p-type 

50 cases. Conclusions are reached in Section 5. 

51 2. Review of Rotating antiferromagnetism theory 

52 2.1. Normal State 

53 We first focus on the normal (non superconducting) state where we review the derivation of rotating 

54 antiferromagnetism (RAF). In section 2.4 we will review the interplay between SC and RAF. Consider 

55 here the t-t' Hubbard model in two dimensions: 



H 



„Cj^a + h.C. 



(1) 
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56 where and j)) are nearest and second-nearest neighboring sites, respectively, and t and t' are 

57 electron hopping energies to nearest and second-nearest neighbors, respectively. Note that hopping to 

58 further neighbors was also considered [5] for more accurate comparison with experiment. The interacting 

59 term in Hamiltonian (1) has been decoupled using 

= (c^Ai) = -(^r) = \Q\e"^\ (2) 

60 and the mean-field theory [1,5-9] was recently combined with the Heisenberg equation [10] in order 

61 to calculate the phase of the this order parameter. To use the Heisenberg equation the interacting term 

62 U rii^nii was rewritten in terms of the spin ladder operators in the following way. In second quantization, 

63 where = c|^Ci_|, the onsite Coulomb repulsion Uui^riii was on one hand written as Urii^riii = 

64 Urii^ — US^S^ and on the other hand as Urii-^Ui^^ = Un-i^^ — US^S^. Summing then dividing by 2 

65 gave the symmetric expression Uni-^riii = ^{rii^ + n^) — ^{S^S^ + S^S^) [10]. The terms S^Sj^ 

66 and Sj^Sl^, which are responsible for onsite spin-flip excitations, contribute by lowering energy for the 

67 sites that are partially occupied by the same density of spin up and down electrons. We decoupled this 

68 term in mean-field theory using (S^) = (c| |Q,i-), which leads to a collective behavior for the spin-flips, 

69 and recovered the results obtained earlier in RAFT [1,5-9]. In this state, a spin flip process at site i 

70 is simultaneously accompanied by another one at another site j; the occurrence of spin flips becomes 

71 synchronized below a transition temperature, which is identified with the pseudogap (PG) temperature. 

72 In Section 2.3 below, an interpretation of RAF in from classical point of view will be given. 

73 The parameter Qi in (2) is thus used to carry on a mean-field decoupling of the t-t' Hubbard model 

74 (1). Consideration of the ansatz where (pi — (pj = n, with i and j labeling any two adjacent lattice sites, 

75 and letting the phase 0j = be site independent but assuming any value in [0, 2n] led to the following 

76 normal state Hamiltonian in RAFT [1,6,7] 

'^i'H'^^. + NUQ^ -NUn\ (3) 

77 where N is the number of sites, and n = {rii^a) is the expectation value of the number operator. Because 

78 of antiferromagnetic correlations a bipartite lattice with sublattices A and B is considered, although no 

79 long-range static antiferromagnetic order is taken into account, and RAFT in only valid away from half- 

80 filling where this long-range order occurs. The summation runs over the reduced (magnetic) Brillouin 

81 zone (RBZ). The Nambu spinor is = (c^| c^|), and the Hamiltonian matrix is 



n 



( -/i' e Qe'^ \ 
e -/i' -Qe^^ 

-Qe~''f' e -fi' J 



E±(k) = V(k)±i?g(k), (4) 



82 yielding the energy spectra 



83 where /i'(k) = n — Un + At' cos cos ky, EqCk) = y e2(k) + (t/Q)^, and e(k) = —2t{cos kx + cos ky). 

84 Using the fact that the energy spectra E± (k) do not depend on phase 0, the matrix "H is transformed 

85 to one that does not depend on using the spin-dependent gauge transformation Cj^i — )■ e'^'^^'^Ci^^ and 
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86 e ^Cj This transformation is equivalent to performing a rotation by angle —0 about the z axis 

87 for the X and y components of the spin operator according to: 



cosd) sin 6 \ I Sf 




sm cos 




Or; 



Note that the thermal averages of S'f and are given by 



i-il = Q COS 0, = -Q sm 0, I e A, or 
n n 

/nx\ /ny\ 

= _gcos0, ^ = gsin0, ieB, (5) 
n h 

and (S'f ) = for i in both sublattices. Because the phase assumes any value between and 2iv, 
rotational symmetry will not look broken for times greater than the period of rotation as we will explain 
below, when we review the calculation of the time dependence of the phase. However if the typical time 
scale of a probe is much smaller than this period symmetry may appear broken. 

The magnitude Q and electron occupation thermal average n are calculated by minimizing the phase- 
independent mean-field free energy. The following mean-field equations were obtained in the normal 
state [1,6,7] 

1 



n 



1 ^ np[E_{k)] - n^E^k)] 

= m • 

96 2.2. Calculation of the time dependence of the phase 

97 The nature of RAF has recently been completely understood after the phase of its order parameter 

98 was calculated as a function of time [10]. Here we summarize how this was done. The Heisenberg 

99 equation = ^ [S^ , H] was calculated in the limit where electron hopping is neglected in comparison 

100 to ^{Sj'S^ + S'Sj'). The values considered in RAFT for onsite Coulomb repulsion are in the range 

101 f/ ~ 3t-5t; this is an intermediate coupling regime where U > t but smaller than the bandwidth ~ 8t 

102 when t' <^ t. Neglecting the effect of electron hopping energies in the Heisenberg equation can be 

103 justified on the ground that spin dynamics is faster than charge dynamics. An onsite spin flip fluctuation 

104 needs a time r ~ to be realized, while a fluctuation caused by a charge hopping between adjacent 

105 sites takes a longer time r ~ (t/ > t). In the Heisenberg equation the bare original interaction was 

106 used instead of RAFT's Hamiltonian (3) in order to treat as best as possible quantum fluctuations. In this 

107 approximation, the following time equation was obtained [10] 

dSt U ^ , 

— ~ i—Sj, T is time, (7) 
rfr Ti ■' 

108 in the intermediate regime where spin dynamics is not governed by the Heisenberg exchange coupling 

109 ~ /U . The latter is suitable for the strong coupling limit. Integration over time r in Eq. (7) gives for 

110 the thermal average 

{Stir)) ^ (5+(0))e^^-/^ (8) 
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111 yielding the phase = Ur/h modulo 2% when (^/(O)) is identified with |(5'+(r))|, (-| (^/(r)) |), for 

112 sublattice A, (B), and e**^ with e*'^^/^. Using this result, the magnetic configuration (5) is rewritten 

113 as follows {Sf)/h = Q cos{usfT), {Sf)/h = ~Q sm{ujsfT) for i in sublattice A or {Sf) /% = 

114 —Qco?,{uJsfT), {Sf)/h = Qsm{usfT) for i in sublattice 5, and (S^) = for i in sublattice A or B. 

115 These thermal averages describe a rotational motion for the spin components with angular frequency 

116 cosf = U /h, and period Tgf = 27rh/U is the time required to perform a spin-flip process, or the time 

117 needed for the rotating order parameter (S^^^^) to complete a 27r-revolution in a classical picture. 

118 2.3. Interpretation of rotating antiferromagnetism 

119 The above derivation of RAF was supported by the following argument that shows that rotating 

120 magnetism (ferro or antiferro) is physically sound and can therefore be realized in a real system 

121 independently of a model. Consider the much simpler case of a single spin precessing in a magnetic 

122 field B along the z-axis, with the initial spin state given by l^'^;, +) = -^{\ t) + I i))- Initially, this 

123 spin points in the positive x-direction. The time-dependent expectation values of this spin's components 

124 are (S'-^) = | cos(a;t), {S^) = | sm{ut), and {S"") = 0, with u = e and mg are the charge and 

125 mass of the electron, and c is the speed of light. This spin's components are therefore confined to rotate 

126 about the z-axis in the xy plane with Larmor angular frequency u. A rotating ferromagnetic state can be 

127 realized by placing such states with the same frequency on a lattice made of N sites. For a rotating 

128 antiferromagetic state, opposite initial states (±|5'j;, +)) where spins point in opposite directions on the 

129 X-axis are placed on any two adjacent sites of a lattice. To relate RAF to spin flip processes, it is noted 

130 that (S^) = (S^) ± i{Sy) = le^*'^* in this example. In a given model, a coupling is necessary for 

131 providing the building bloc for RAF, which is a spin precessing about an effective magnetic field (with 

132 no local magnetization) for each lattice site and the anti-alignment of the adjacent rotating moments. The 

133 RAF state constructed in this way shows a hidden order that can be realized even at finite temperature 

134 without violating the Mermin-Wagner theorem [11]. 

135 The above simple case allowed us to interpret RAF as a state where spins precess collectively in a 

136 synchronized way in the spins' xy plane around a staggered effective magnetic field B = mecU/h\e\ 

137 generated by onsite Coulomb repulsion, h/2 in (S*^) = |e^*'^* is replaced by the magnitude of the RAF 

138 order parameter Qh, which assumes values smaller than h/2 due to thermal fluctuations. In comparison 

139 to ordinary spin waves in an antiferromagnet, RAF's state was interpreted as a single q = (vr, tt) spin 

140 wave occuring as a consequence of zero staggered static magnetization. The spin-wave theory does 

141 not however apply for our system (where (S*/) = 0), since this theory has to be built on top of a Neel 

142 background with finite (S^). Also, in comparison to ordinary antiferromagnetic spin-density order, RAF 

143 is characterized by a local magnetization that is not static because of the time dependence of the phase 

144 of the magnetization. It is thus clear that RAF will have all the typical effects of spin-density order on 

145 the evolution of the electronic structure with doping, but is expected to go undetected for experimental 

146 probes like neutrons due to the time dependence of the phase. We predicted [10] that rotational symmetry 

147 will not look broken for experimental probes that are characterized by a time scale greater than the period 

148 of rotation Tgj = 2Txh/U of the rotating order parameter of RAF. For such probes, averaging over times 

149 longer than the period will not allow for the observation of RAF. In RAFT, electron hopping enery t 
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150 taken to be 0.1 eV in fitting data. Taking U = 3t = 0.3 eV gives Tgf ^ 10 s. For neutrons for 

151 example the typical time would be the time spent by a given neutron in the immediate vicinity of a given 

152 spin during the scattering process. If this tme is greater than the period Tgf then neutrons will not detect 

153 RAF. If the time spent by the neutron in the vicinity of the spin is smaller then there is a chance RAF 

154 will be detected. Note that smaller times means higher energies for neutrons. This is issue is still under 

155 investigation and will be reported on when this investigation is over. 

156 2.4. Interplay between RAF and SC 

157 In RAFT, d-wave SC was introduced phenomenologically using an attractive coupling between 

158 electrons on adjacent sites. The term — '^i.t'^j,^ 1^°^ added to Hamiltonian (1), and is 

159 decoupled using -D(ij) = (Q^;Cj^-|-). To get a d-wave gap we set -D(ij) = Dq along the x-direction 

160 and -D(jj) = —-Do along the y-direction [1,6]. 

161 When both SC and RAF orders are taken into account, the mean-field Hamiltonian is written in terms 

162 of an eight-component spinor given by 

163 and assumes the expression [1,6] 

H = ^k^^k + UNQ"^ + UNm'^ 

k< 

+AVNDI - UNn^ - fiN, (10) 

164 where "H is an 8 x 8 matrix: 

H = 

165 with H' and Uq, two 4x4 matrices, given by 



166 and 

/ gf/ \ 

-QU 
^ -QU 

^ QU J 

167 The k-dependent superconducting gap is D(k) = 2yDo(cos — cos ky) where Dq = (c^ij,4.c^i±iji-) = 

168 — (c^2j4,c^2j_ti,t) involves two adjacent sites on different sublattices. Triplet SC [12] is ruled out by 

169 choosing (c^i j,4.cfj_,_]^ ^l) = — (c^ij-|-cfj_,_]^j l); so only spin-singlet SC, which is relevant to HTSCs, is 

170 considered in RAFT. The way decoupling is done using creation and annihilation operators rather than 

171 the spin- singlet and triplet superconducting operators, which are a combination of products of the c's 

172 and c^'s allowed us to avoid the generation of tripet SC if it is not present initially [12]. 



-Uq -w ) 

e(k) D(k) \ 
-fi'{k) D{k) 

D{k) ^'(k) -6(k) 

-e(k) /i'(k) i 
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173 In Hamiltonian (10), e(k) = —2t{coskx + cos ky), and /i'(k) = ji — Un + 4t' cos /c^; cos /c^ have 

174 the same expressions as in the absence of SC. In equation (10), the summation over k takes into 

175 account the doubling of the Brillouin zone and the fact that summation is now over k and — k. The 

176 size of the mean-field Hamiltonian matrix 1-L is twice as large as that in the density rf-wave (DDW) 

177 approach [13], proposed for the PG behavior, or in other approaches that deal with the interplay between 

178 antiferromomagnetism and SC [12]. 

179 The energy spectra obtained by diagonalizing the matrix "H are ±£'i(k) and ±_E'2(k) with 



£;,(k) = J[/i'(k) + (-i)-E,(k)p + D2(k), ,. = 1,2, (11) 



180 where ^<,(k) = y'e2(k) + QW^. 

181 Minimizing the free energy function with respect to Q and Dq, and calculating the density of electrons 

182 n led to the following mean-field equations that describe the interplay between RAF and SC for HTSCs 

183 with tetragonal symmetry: 



n = -^ E A.tanh(^) + -, (12) 

k, 1^=1, 2 

184 where 

AM = [-/^'(k) - {-irEM]/EM- (13) 

185 In the case of crystals with orthorhombic symmetry it is possible that Di^ij)^ ^ ~E>{i,j)y because the 

186 superconducting coupling constants along the x axis and Vy along the y axis may differ. Then, 

187 the superconducting gap takes on the form D(k) = tps(coskx + cos ky) + ipd^cosk^ — cos ky) with 

188 ips = VxDx — VyDy, and ipd = V^D^ + VyDy, which implies that it shows rf-i-s-pairing symmetry. This 

189 is a consequence of the absence of invariance under rotations by 7r/2 of the CUO2 plane, and is therefore 

190 consistent with arguments based on group theory [14]. 



191 3. RAF and SC parameters versus temperature and doping 

192 The mean-field equations (12) are solved numerically in order to get the parameters Q, Dq, and 

193 average number of electrons per site and per spin n. In hole-doped (p-type) systems, the density of holes 

194 given hy p = 1 — 2n is the density of electrons missing below half-filling. For electron-doped (n-type) 

195 systems the density of electrons rig = 2n — 1 is the density of electrons above half-filling. Note that at 

196 half-filling n = 1/2, so that Ue = p = in this case. 

197 Fig. 1 shows RAF and SC parameters Q and Dq versus doping for the Hamiltonian parameters 

198 U = 2.8t, V = 0.85t, and t' = — 0.16t. If RAF were not taken into account, SC would be optimum 

199 at half filling. When RAF is allowed in, SC is destroyed near half filling, but coexists with RAF in the 

200 underdoped regime for p-type doping, and for all doping values where -Dq 7^ in n-type case. The doping 
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Figure 1. The typical behavior of the rotating order parameter Q and superconducting 
parameter Dq with doping is illustrated here for U = 2.8t, V = 0.85t, and t' = — 0.16t. 
Temperature is T = 0.05t. The behavior shown is practically the same at zero temeprature. 
The doping values where Q vanishes in both p and n-type cases are interpreted as QCP. 




Figure 2. The temperature dependence of Q and Dq is displayed for two values of doping in 
the p-type case. The Hamiltonian parameters are U = 2.8t, V = 0.85t and t' = — 0.16t. 
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Figure 3. The temperature dependence of Q and -Do is displayed for an n-type case with 
Ue = 0.075. The Hamiltonian parameters used are f/ = 2.8t, V = 0.85t and t' = — 0.16t. 
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201 values where Q vanishes in both p and n-type systems are identified as quantum critical points (QCP) [ 1 ] . 

202 In Fig. 1 , the QCP occurs within the superconducting dome in the p-type system, but outside of the dome 

203 for the n-type system. RAF's parameter Q has been proposed to model the PG in HTSCs [1,6], and the 

204 PG temperature T* has been identified with the temperature below which Q becomes nonzero. Figures 

205 2 and 3 display the temperature dependence for Q and Dq for some given doping levels. Using these 

206 kind of figures, the PG (T*) and SC (Tc) temperatures were calculated in Ref. [6] for p-type cuprates. 

207 Again the competition is apparent between SC and RAf in the p-type case, because as soon as Do 

208 becomes nonzero RAF's order parameter Q decreases monotonically as seen in Fig. 2 for p = 0.1 in the 

209 underdoped regime. Note that the optimal doping for the Hamiltonian parameters used here is p ~ 0.20 

210 in the p-type case, and ~ 0.075 in the n-type case; see Fig. 1. For p = 0.195 (close to optimal 

211 doping), it is interesting to note that Q decreases significantly below Tc when Dq becomes nonzero. The 

212 behavior for n-type case is totally different. For Ue = 0.075 (optimal doping), Q barely decreases when 

213 Dq becomes nonzero below Tc, then even increases slightly and saturates at low temperature as the inset 

214 of Fig. 3 shows. The phase diagram obtained by letting in a Bardeen-Cooper-Schrieffer [15] picture 

215 T* ~ Q(T = 0) and Tc ~ Do{T = 0) in Fig. 1 is in qualitative agreement with experiment for p-type 

216 La2_KSra;.Cu04 and n-type Nd2-xCea:Cu04 cuprates. For the latter, the PG is reported to vanish almost 

217 when SC does in the overdoped regime [2]. For the former significant experimental evidence suggests 

218 the disappearance of the PG rather closer to the optimal doping [16]. 

219 4. Doping dependence of electronic structure 

220 4.1. Analysis of energy spectra 

221 As mentioned in the previous section, the appearance of RAF below a critical value of doping as the 

222 latter is reduced from overdoped to underdoped regime for p-type or n-type systems at zero temperature 

223 has been interpreted as a QCP. The case of p-type has been discussed before [1,6,8]. This QCP induces a 

224 reconstruction of the FS practically in the same way an ordinary spin-density order does [17]. However, 

225 RAF is not an ordinary spin density order as explained in Section 2. Figure 4 shows energy spectra along 

226 symmetry lines of the Brillouin zone for T = O.lt. For p = 0.075, one clearly sees a gap at (vr, 0), in 

227 agreement with experimental data for La2-a;Sr3.Cu04 [18]. Also a small hole-like band is seen along the 
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Figure 4. The energies E^(k) and -E+(k) are plotted versus k along symmetry lines of the 
Brillouin zone. The Hamiltonian parameters are f/ = 2.8t, V = 0.85t and t' = — 0.16t, and 
hole doping is p = 0.075 in the underdoped phase for the figure on the left, and p = 0.24 
in the ovedoped regime for the figure on the right. The dotted horizontal line indicates the 
position of the Fermi energy. Temperature is T = O.lt and Dq = 0. 




0.5 0.5 0.5 1 0.5 0.5 0.5 1 

(0,0) (7t,0) (7t,Jl) (0,0) (0,0) (71,0) (71,71) (0,0) 



Figure 5. The energies -E'-(k) and -E'+(k) are plotted versus k along symmetry lines of the 
Brillouin zone. The Hamiltonian parameters are f/ = 2.8t, V = 0.85t and t' = — 0.16t, and 
electron doping is = 0.06 in the underdoped phase for the figure on the left, and = 0.2 
in the ovedoped regime for the figure on the right. Temperature is T = O.lt and -Do = 0. 




""•■O 0.5 0.5 0.5 1 "^0 0.5 0.5 0.5 1 
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228 diagonal around (7r/2, 7r/2). The origin of the gap at (tt, 0) for this doping and the small hoe-like band 

229 in the vicinity of (vr/2, n/2) is the nonzero value of RAF's order parameter Q. This gap is responsile 

230 for the PG behavior in the underdoped regime. The hole-like band is also seen along the RBZ bounday 

231 [(tt, 0) (0, 7r)] as shown in Fig. 6 for p = 0.1. For p = 0.24 in the overdoped regime, the PG has 

232 closed and the hole-like pocket has reached the (tt, 0) and (0, vr) points as can be seen along the RBZ 

233 boundary in Fig. 6. Along this boundary -E'+(k) = E_(k), when Q = at T = O.lt, is above the 

234 chemical potential all the way between (vr, 0) and (0, tt) . 

235 For the n-type case. Figure 5 displays the spectra for doping Ue = 0.06 in the underdoped regime, and 

236 for rig = 0.2 well in the overdoped regime for T = 0. It. The PG behavior is now a consequence of a gap 

237 at (7r/2, 7r/2), and a small electron pocket forms near (tt, 0). For Ue = 0.2, the PG is zero because Q has 

238 vanished, and the electron pocket at (vr, 0) joined that at (0, tt) . This can be understood by examining the 

239 spectrum along the RBZ boundary which gives a comletely full band along this direction. For example, 

240 for rig = 0.21 in Fig. 6, E+(k) = E^(k) < 0, which means that these bands are full. The above analysis 

241 can be made even more transparent by calculating the FS, a task undertaken below. 
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Figure 6. The energy spectra along the boundary of the magnetic Brillouin zone ((tt, 0) — >■ 
(0, tt)) is shown. Temperature is T = O.lt and Dq = 0. 




242 4.2. Evolution of the Fermi surface with doping 

243 In RAFT, the occupation probability n(k) was defined by writing the average number of electron per 

244 spin and site n in equation (12) as n = X^k '^(k), which yields [1] 

n{k) = -\j: am tanh[^M^] + 1 (14) 

245 A^(k) is given in Eq. (13). n{V) was then interpreted as the probability that the state with wave vector k 

246 is occupied by an electron with spin up or down. 

247 Ronning et al. [19] extracted n(k) by integrating ARPES energy distribution curves over energy for 

248 the material Ca2Cu02Cl2, then deduced the FS by locating the steepest drops in n(k) in analogy with 

249 a Fermi gaz. Also, using the same method the Fermi surface for Bi2Sr2CaCu208+5 in the overdoped 

250 regime was obtained. Here we implement the same argument in RAFT, namely the FS is determined by 

251 the sharp drops in the occupation probability. This method was also used in Ref. [5] and gave results in 

252 agreement with the determination of the FS using the spectral function. 

253 Fig. 7 shows two-dimensional plots of n(k) for three doping levels in the p-type case for Hamiltonian 

254 parameters \J = 2.8t, V = 0.85t, and t' = — 0.16t at temperature T = O.lt. This temperature is 

255 above any transition temperature for SC. The FS is made of hole pockets around (±7r/2, ±7r/2) in the 

256 underdoped regime as shown for p = 0.06. The energy spectra in Figs. 4 and 6 show well that in the 

257 presence of the PG, the upper E_^. and lower _E„ bands are separted by gaps along all the symmetry 

258 lines in the underdoped regime. Hole-like pockets can clearly be seen for p = 0.075 in Fig. 4 around 

259 (vr/2, 7r/2). Around optimal doping p = 0.2, the hole pockets reach the points (ivr, 0) and (0, ±7r). In 

260 the ovedoped regime, where the PG is zero, the FS is made of large contours around (0, 0) and (tt, tt) 

261 as can be seen in Fig. 4 for p = 0.24. For the latter, because the PG is zero the upper band and 

262 lower band touch at (±7r, 0), (0, ±7r) and (±7r/2, 7r/2) to form a tight-binding spectrum given by 

263 -E'±(k) = ±\2t{cos kx + cos ky)\ — At' cos cos ky — fi + Un. The presence of the absolute value in this 

264 tight-binding energy is a consequence in the limit Q — )■ in y^e^(k) + U'^Q'^ for the overdoped regime. 

265 In the p-type case, the FS in RAFT thus evolves strongly with doping. It reconstructs at the QCP doping 

266 where Q vanishes, its topology changes from small hole-like pockets in the underdoped regime below 

267 this QCP to large contours in the overdoped regime. This is qualitatively consistent with the quantum 

268 oscillations observed in resistivity by Doiron-Leyraud et al. [20] for YBa2Cu306.5, which indicated that 



Version December 12, 2012 submitted to Symmetry 



12 of 15 



Figure 7. The occupation probability ?7,(k) is shown in the Brillouin zone. The hole densities 
are shown on the graphs. The Hamiltonian parameters are U = 2.8t, V = 0.85t, and 
t' = -O.m, andT = O.lt. 




269 a well defined small FS characterizes this underdoped cuprate. Subsequent work by Sebastian et al. [21] 

270 for YBa2Cu306+x supported the existence of small closed pockets in the underdoped regime as well. 

271 The calculated FS undergoes also a significant reconstruction when doping changes from p-type to 

272 n-type across half-filling, (note that RAFT is only valid outside of the HTSCs' AF phase around half- 

273 filling). In RAFT, for the Hamiltonian parameters considered here the FS in the underdoped regime 

274 for n-type cuprates consists of electron pockets around points (±7r, 0) and (0, ivr), rather than pockets 

275 around (±7r/2, ±7r/2) in the underdoped regime of p-type cuprates. This is clearly seen in Fig. 8 for 

276 Tie = 0.06 and = 0.1, and is consistent with the energy spectra in Fig. 5, which show well the 

277 existence of a small electron pocket at (tt, 0) for rie = 0.06. When the PG vanishes in the overdoped 

278 regime, the electron pockets join to form large contours as seen for = 0.2. Armitage et al. [3] reported 

279 ARPES data for n-type material Nd2-a;Ce2;Cu04 which can be interpreted as revealing the existence of 

280 pockets around (tt, 0) and symmetric points in the underdoped regime. Also, Matsui et al. [22] measured 

281 the evolution of the FS with doping for this material using ARPES. A close look at Fig. 1 of their work 

282 reveals a FS mainly near k-points (tt, 0) and (0, tt) for doping x = 0.13, but the FS evolves into larger 

283 countours joining these two points for the larger doping levels x = 0.16 and x = 0.17. Note that if 

284 one symmetrizes Matsui et a/.'s FSs about the line joining (tt, 0) and (0, tt), one will get FSs that look 

285 similar to those calculated here, and shown in Fig. 8. Their measurements were done only along the 

286 FS in the thight-binding limit. We predict that if measurements were performed along the image of this 

287 tight-binding FS with respect to line (tt, 0)-(0, tt), then one would obtain a FS that looks like ours. Also, 

288 the reconstruction of the FS, as illustrated in Fig. 8, at the QCP doping where Q vanishes (so where the 

289 PG vanishes) is in agreement with SdH oscillation results of Helm et al. [2] for the above material. SdH 

290 oscillations revealed a FS evolving from small pockets to large contours as doping goes from underdoped 

291 regime to overdoped regime. 
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Figure 8. The occupation probability n(k) is shown in the Brillouin zone. The electron 
densities above half-filling Ue are shown on the graphs. The Hamiltonian parameters are 
U = 2.8t, V = 0.85t and t' = -0.16t and temperature is T = O.lt. 




292 5. Conclusions 

293 We studied the reconstruction of the Fermi surface under the effect of hidden rotating antiferro- 

294 magnetic order in both p-type and n-type high-Tc; cuprates. For the Hamiltonian parameters used 

295 here, as a consequence of rotating antiferromagnetism, the Fermi surface in p-type case reconstructs 

296 at the quantum critical point near optimal doping where the pseudogap vanishes. It consists of hole 

297 pockets around (±7r/2, ±7r/2) in the underdoped regime, but changes to large closed contours in the 

298 overdoped regime. For n-type cuprates, the location and topology of the Fermi surface is different than 

299 in p-type materials. The hole pockets are of a form resembling squares around (±7r, 0) and (0, ±7r) 

300 in deep underdoped n-type systems. When the pseudogap becomes zero beyond the quantum critical 

301 point in the ovedoped regime, the Fermi surface changes to large closed contours. These results are in 

302 good qualitative agreement with experimental data for La2-a;Sr3;Cu04 and Nd2-xCea;Cu04 for the set of 

303 Hamiltonian parameters used in the present calculations. 
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